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Abstract: We discuss how to consistently use Effective Field Theories (EFTs) to set uni¬ 
versal bounds on heavy-mediator Dark Matter at colliders, without prejudice on the model 
underlying a given effective interaction. We illustrate the method for a Majorana fermion, 
universally coupled to the Standard Model quarks via a dimension-6 axial-axial four-fermion 
operator. We recast the ATLAS mono-jet analysis and show that a considerable fraction 
of the parameter space, seemingly excluded by a naive EFT interpretation, is actually still 
unexplored. Consistently set EFT limits can be reinterpreted in any specific underlying 
model. We provide two explicit examples for the chosen operator and compare the reach 
of our model-independent method with that obtainable by dedicated analyses. 
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1 Introduction 

After the discovery [1, 2] of a Higgs boson compatible [3, 4] with the Standard Model (SM), 
and the non-detection so far of new particles [5] at the LHC, searches for Dark Matter 
(DM) in the form of a Weakly Interacting Massive Particle (WIMP) are becoming a central 
theme for the LHC general purpose experiments (for a pedagogical review, see e.g. [6]). 

In the WIMP hypothesis, the DM particle has a mass in the GeV-TeV range, and 
the strength of its couplings to the SM particles is roughly of electroweak size. The relic 
density, generated by thermal freeze-out, can then match the cosmological and astrophysical 
observations. This so-called “WIMP miracle” receives further support from the fact that 
WIMPs are ubiquitous in new physics models of Electro Weak Symmetry Breaking (EWSB), 
motivated by the naturalness problem of the SM. The latter is an appealing and intensively 
explored possibility, but WIMP DM might well originate from a completely unrelated sector. 
Moreover, we currently have no idea of how the complete EWSB sector looks like, thus there 
is not much we can say a priori on the specific properties of WIMP DM. 

In the situation described above, a general and model-independent exploration appears 
mandatory. Commitment to specific benchmark models (or classes of benchmark models) 
should be avoided whenever possible in the analysis of experimental data, or at least treated 
as an accessory step in the interpretation. The goal is to search for WIMP DM in a 
comprehensive way, leaving no unexplored corners in theory space. ^ 

^For example, when planning future direct-detection experiments sensitive to low-mass WIMPs in the 
1-10 GeV range, it may be important to understand on general grounds how much room for discovery is 
left after the so far unsuccessful LHC searches, without being committed to specific benchmark models. 
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In the case of heavy-mediator DM, this program can be carried out, at least to some 
extent. The working hypothesis is that the DM candidate X interacts with the SM through 
the exchange of one or more particles, called “mediators”, whose mass is well above the mass 
mDM of the DM particle. This assumption is motivated by the present lack of evidence for 
new particles at the LHC, but it is not the only possibility. The case in which the mediator 
is a SM particle, such as a weak or the Higgs boson, is equally plausible and deserves equal 
attention. Light and very weakly coupled mediators can be also conceived. 

Focusing on the heavy-mediator case for the rest of this paper, it is relatively easy to 
set up a model-independent strategy for DM searches. We can exploit the fact that the 
dynamics of the DM particle can be universally described, in the appropriate kinematical 
regime, by a low-energy EFT Lagrangian [7-20], invariant under the SM gauge group and 
the Lorentz group:^ 

£eFT = >CsM + + Ant • (1-1) 

In the above equation, Tsm denotes the SM Lagrangian, Cx is the free Lagrangian for X, 
and Tint contains the operators describing the DM interactions with the SM particles, plus 
possible additional interactions in the DM and SM sectors. If we knew the true microscopic 
DM theory, these operators could be computed by integrating out the mediators. However, 
their form is universal and the lack of information on the mediator dynamics merely prevents 
us from computing the value of their coefficients, which are thus free input parameters of 
the EFT. 

The allowed operators in Tint can be classified according to their mass dimension d, 
for different hypotheses on the DM quantum numbers. In many relevant cases the DM 
quantum numbers forbid renormalizable interactions with d < 4, and the lowest-dimensional 
operators have d = 5,6. For the physics to be considered in this paper, we can assume that 
the d = 5 operators are negligible and the leading operators have d = 6: 


Tint — 


1 




( 1 . 2 ) 


where the sum runs over all d = 6 operators Oi allowed by the symmetries, a are di¬ 
mensionless coefficients and the overall effective coupling strength is parameterized by a 
dimensionful coupling 1/M^. 

While the EFT can be formally defined independently of any consideration about its 
microscopic origin, its range of applicability and thus its physical relevance depend on the 
underlying theory. Namely, the EFT provides an accurate description of the underlying 
model only for elementary scattering processes taking place at a low enough centre-of-mass 
energy Ecm, below a certain critical scale Mcut usually called the EFT cutoff. This cutoff is 
determined by the mass of the mediators in the microscopic theory but it is unknown from 
the viewpoint of the EFT and it should thus be treated as a free parameter, on the same 
footing as those introduced above. 


^At energies as low as those relevant for direct detection experiments, it may even be convenient to switch 
to a non-relativistic EFT [21-24], but for obvious reasons this approach precludes a direct comparison with 
collider searches and will not be pursued here. 
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The EFT is then characterised by at least three parameters: 

• the DM mass mnivi; 

• the scale M* of the interaction; 

• the cutoff scale Mcut- 

If a single operator appears in eq. ( 1 . 2 ), the corresponding dimensionless coefficient can be 
absorbed in M*, otherwise the EFT parameters also include the Cj coefficients. With these 
free parameters, the EFT faithfully reproduces the predictions of any microscopic theory 
for all processes taking place at i?cm < Mcnt- Given that the effective operators in eq. (1.2) 
may have many possible microscopic origins, exemplified by the plethora of models in the 
literature, this simplification is particularly useful. 

Notice that Mcut and M* are logically independent parameters, however they can be 
approximately related by 

Mcut = g*M^, (1.3) 

where 5 * is a suitably defined coupling strength of the underlying microscopic theory. The 
simplest way to motivate the above equation is the analogy with the Fermi theory of weak 
interactions, where the cutoff Mcut is the mass of the W boson (the mediator in this 
context), gf* is the SU{2) gauge coupling and 1/M^ is the Fermi constant Gp: they 
indeed obey eq. (1.3) up to numerical factors. Alternatively, the physical meaning of g^ can 
be appreciated by noticing that the EFT interaction strength is given, for processes taking 
place at a given center-of-mass energy Ecm, by the dimensionless combination 
At the mediator scale, i.e. the cutoff scale Mcut, this strength becomes M^^^/M^ = g'^, 
providing further justification for interpreting g^ as the typical mediator coupling. Using 
eq. (1.3) to re-express M* in terms of 5 * will be important in section 2.2, in order to draw 
the current limits on a plane suited for theoretical interpretation. 

The EFT can be straightforwardly used to predict the cross-sections for a number of 
relevant reactions, namely the DM annihilation in the Early Universe, which determines 
the thermal relic density, the present-day annihilation, which controls indirect detection, 
and the DM scattering on nucleons, which direct search experiments try to detect. Indeed, 
all these reactions take place at safely small Ecm and therefore, up to subtle effects that 
might be encountered in the relic density calculation, the EFT predictions are automati¬ 
cally trustable. If collider searches could be added to the list, we would reach the truly 
remarkable conclusion that all the experimental information on heavy-mediator DM can be 
simultaneously interpreted and compared in a completely model-independent fashion, with 
no prejudice on the specific nature of the mediator and of its couplings to DM and to the 
SM. However, the usage of the EFT at colliders is problematic, because the energy of the 
reaction in which the DM is produced is not necessarily smaller than Mcut, and this risks to 
invalidate the EFT predictions. The effect is quantitatively amplified by the requirement 
of extra hard objects (e.g., one jet), in addition to the pairs of DM particles, for the signal 
to be triggered and disentangled from the background. This problem has been discussed 
at length in the recent literature (see e.g. refs. [25-34]), the goal of the present article is to 
illustrate a simple and practical solution. 
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The basic observation is that the processes for DM production at colliders can be split 
into two kinematically distinct classes, characterised by a centre-of-mass energy below and 
above Mcut, respectively. The former class defines our theoretical signal, and its rate is 
accurately predicted by the EFT. The latter would instead require the knowledge of the 
microscopic theory and its contribution to the cross-section is thus unpredictable within 
the EFT. Under certain conditions, to be described below, the second class can be simply 
ignored and an experimental limit can be set on the signal defined, as explained above, by 
the DM production reaction restricted to £'cm < Mcut- This is possible if the experimental 
search is performed as a counting experiment in one or several signal regions, defined by 
a certain set of cuts on the visible final state particles. The low and high Ecm processes 
both contribute to each signal region, but in a purely additive way, since low and high Ecm 
regions are quantum-mechanically distinguishable and do not interfere. Therefore a lower 
bound on the expected cross-section is obtained by considering only the “well-predicted” 
signal events, namely those restricted to the Ecm < Mcut region. If the result of the search 
is negative, an exclusion upper bound fiexc is set on the cross-section, which we can interpret 
through the inequality 


s 

cteft 




— '^true CTexc , 


(1.4) 


where denotes the “true” signal as it would be computed in the unknown microscopic 
theory. Regardless of what the latter theory is, the restricted EFT signal cteft system¬ 
atically underestimates the cross-section and thus provides a conservative, but correct, 
exclusion limit. ^ 

The rest of the paper is organised as follows. In section 2 we illustrate our limit-setting 
strategy in the explicit example of a four-fermion operator obtained as the product of axial 
currents involving the SM quarks and a SM-singlet Majorana fermion DM. This choice is 
partly motivated by the fact that direct and indirect detection experiments have a poor 
sensitivity to this operator, thus collider searches are expected to be the most sensitive 
ones, but the same method can be applied to all other operators. We quantify the reach of 
current collider searches by recasting the ATLAS mono-jet results available at the time of 
this work, and show how the latter can be presented in a theoretically useful way. Besides 
the methodological proposal, the important physics point is that, from the general EFT 
viewpoint, the present collider bounds on DM have not yet probed the most plausible region 
of parameter space. To access such region, we need not only more energy and luminosity, as 
expected in the forthcoming runs of the LHC, but also improvements in the experimental 
analyses. In section 3 we describe another relevant feature of our strategy, the fact that the 
limits set in the EFT can be straightforwardly re-interpreted as constraints on any specific 
microscopic model. This is because the EFT parameters can be computed in the underlying 
microscopic theory and expressed in terms of the “fundamental” parameters of the latter (for 
previous discussions of the interplay between EFT and underlying microscopic models in 
DM searches at colliders, see again refs. [25-34]). We consider two representative models. 
Model A and Model B, which both give rise to the same axial-axial effective operator, 
and compare the limits derived from the EFT with those obtainable from a dedicated 


®For a similar approach in the context of Higgs EFTs, see [35]. 
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interpretation of the mono-jet search within the two models. Since our signal cross-section 
systematically underestimates that of the microscopic theory, we obtain conservative limits. 
We find that these limits differ significantly from those obtained in the full models only in 
the kinematical region where the mediators can be resonantly produced. In such a case, 
however, more comprehensive experimental strategies, complementing the event selection 
used for heavy-mediator DM searches with other selections that can take full advantage 
of the resonant production of the mediators (single or in pairs, with one or more jets in 
the event), should be able to provide stronger bounds. We end this section by discussing 
two aspects of our simple and practical approach that can be helpful for the comparison 
with a similar but more model-dependent approach recently put forward in [27, 29, 30]. We 
finally present our conclusions in Section 4. Some back-up material is collected in three 
appendices. Appendices A and B provide details on Model A and Model B, respectively. 
Appendix C collects the approximate analytical formulae used to draw the relic density 
constraint in some of the figures. 

2 Limit-setting strategy 

For the present study, we assume that the DM particle is a Majorana fermion, singlet under 
the SM gauge group and represented by a self-conjugate four-component spinor X = X^, 
whose free Lagrangian reads 

Cx =-X {i(^ — rriDu) X . (2.1) 

As for the interactions between X and the SM particles, we just choose a representative 
example to illustrate our limit-setting strategy, assuming that they can be described, in the 
low-energy limit, by the single ^ axial-axial four-fermion operator ^ 

, (2 

where the sum is over all quark flavours {q = u, d, c, s, t, b), the dimensionless coefficient c has 

been re-absorbed in the definition of M*, and the overall minus sign is purely conventional 

in the present context. This effective operator mediates DM pair-production at the LHC, 
a process which is however undetectable and impossible to trigger because of the lack 
of visible objects in the final state. Searches are performed by considering extra visible 
emissions from the initial quarks, leading to the so-called “mono-D” signatures, where V 
could be a jet [38-43], a photon [44-47], a massive weak boson [48, 49] or a top quark 
[50, 51]. Below we restrict our attention to the mono-jet searches, because they currently 
show the best sensitivity, but our considerations also apply to the other channels. 

Radiative corrections may generate additional operators [36, 37], this can be important when comparing 
with direct dark matter searches but does not play a role in the present context. 

®This operator is twice the M6 operator in [12], and formally coincides with the D8 operator in [15], 
which is often taken as a benchmark for experimental searches. Notice however that we are dealing with a 
Majorana spinor normalised as in (2.1), while D8 involves a canonically normalised Dirac spinor. 
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signal region 

SRI 

SR2 

SR3 

SR4 

and .B™'®® 

>120 

>220 

>350 

>500 

Texc[pb] 

2.7 

0.15 

4.810-2 

1.510-2 


Table 1. Signal region definitions (cuts expressed in GeV) and 95% CL limits from ref. [42]. 

2.1 ATLAS mono-jet recast 

Searches for a jet plus missing transverse energy (if™'®®) have been performed at the LHC 
by the ATLAS [39, 41, 42] and CMS [38, 40, 43] collaborations.® We focus here on the 
ATLAS analysis in ref. [42] because, among those available at the time of the present work, 
it is particularly suited to illustrate the general point we would like to make. The search 
is performed as a counting experiment in four overlapping signal regions (SR), with pre¬ 
selected events characterized by if™'®® > 120 GeV, one jet with > 120 GeV, jr^j < 2 and 
at most one additional jet withpx > 30 GeV and \r]\ < 4.5. If found, the second jet is asked 
to be separated in the azimuthal direction from the pip"®® by a cut A(/> >0.5. Additional 
requirements, namely on the primary vertex reconstruction and on the absence of extra jets 
with anomalous charged/calorimetric composition, are not directly relevant for our study, 
since their impact crucially depends on the detector response, which we cannot simulate. 
The four signal regions SRi {i = 1, 2, 3,4) are defined by increasingly strong cuts on if™'®® 
and on p!^^. The results are presented as upper bounds, on the visible cross-section in 
each region. The SR definitions and the exclusion limits are summarized in table 1. 

We reinterpret these limits as follows. The expected signal in each SR is expressed as 

o-sRi = a X AiX et, (2.3) 

where a denotes the total signal cross-section defined as in eq. (1.4), Aj is the geometric 
cut acceptance, as obtained from a leading-order parton-level simulation, and the efficiency 
€i is the correction due to showering, hadronization and detector effects. Acceptances 
and efficiencies depend on the DM mass tudm and on the cutoff Mcnt, while the operator 
scale M* only enters the total cross-section as an overall factor 1/Mf. We compute the 
parton-level quantities a and Aj by MadGraph 5 [53] simulations, while we estimate the e* 
corrections by matching with the limits on the D 8 operator scale reported in ref. [42]. In 
practice, we simulate the same D 8 operator signal considered in ref. [42] (i.e. Mcnt = oo 
in eq. (1.4)), we compute a x Ai and we determine e* such as to reproduce the ATLAS 
limit on the effective operator scale as a function of the DM mass. Actually, since only the 
third SR is used by ATLAS to set the limit, only 63 can be obtained in this way. The same 
efficiencies are used for the other SRs, although we see no reason why the efficiency should 
stay the same in all the regions. The result of this procedure gives rather small efficiencies, 
of around 60%, approximately constant over the whole DM mass range. We verified that 
this considerable signal loss is mainly due to the fact that our simulation does not include 
the showering-level production of extra jets, a significant fraction of which are vetoed in 

®For a very recent update of the ATLAS mono-jet analysis, which appeared after the completion of our 
work, see also [52]. 
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the ATLAS event selection. Notice that the efficiencies for our signal might be significantly 
different from those estimated in the naive EFT because, although based on the same 
effective operator D8 of eq. (2.2), our signal is constrained by Mcut to the low invariant 
mass region, thus it is expected to have different kinematical distributions. A complete 
simulation in different regions of ut-dm ^-nd Mcut, including showering and matching, would 
be needed for an accurate analysis, but goes beyond the aim of the present illustrative 
example. 

Under the assumptions explained above, the expected signal takes the form 

o'SRi(M*,mDM,Mcut) = cr(M*,mDM,Mcut) X Ai(mDM,Mcut) X e 

1 T'gV^ " 

= X CT(mDM,Mcut) X Aj(mDM,Mcut) X e, (2.4) 

where the overall scaling of the cross section with M* has been factored out and the result 
expressed in terms of a reference cross-section a computed for M* = 1 TeV. The reference 
cross-section times the acceptances are obtained by MadGraph 5 [53] simulations of DM pair 
plus one parton production, duly restricted by the hard jet kinematical cuts that define 
each SR. El™'®® cuts are automatically imposed because the jet and the missing transverse 
momentum, i.e. the transverse momentum of the DM pair, are back-to-back in our parton- 
level sample. The theoretical restriction Ecm < Mcut, which ensures the validity of the 
EFT description as explained in the introduction, should be imposed as a cut on the total 
invariant mass of the hard final states of the reaction, namely as 

[p(DMi) + p(DM2) + pi''*]' < . (2.5) 

For our parton level simulation this is equivalent to a cut \/s < Mcut on the total partonic 
centre-of-mass energy, however when going to the showered and matched level one should 
be careful not to cut on \/^ but on the variable in eq. (2.5), with pi®* the leading jet 
four-momentum. 

A scan is performed in the (mDM,Mcut) plane for each SR and the values of cr x Aj 
are used to construct two-dimensional interpolating functions. A significant dependence 
on ttidm is only found for tudm ^ 80 GeV, while for smaller values a x Ai is basically 
constant in tudm- Once the signal cross-sections are known, the 95% CL limits are imposed 
as constraints 

'7SRi(M*, ruDM, Mcut) < o-^c ) (2-6) 

out of which the 95% CL allowed regions are immediately found in the three-dimensional 
parameter space (M*, tudm, Mcut)- The limits from the various signal regions can be studied 
separately or combined. For our illustrative purposes, the combination will be performed by 
just taking the overlap of the four allowed regions. The results of this simple limit-setting 
procedure are discussed in the following section. 

2.2 Results and discussion 

At fixed moM and Mcut, the ATLAS limits in eq. (2.6) become lower bounds on the scale 
M*, reported in fig. 1 as functions of ttidm and for different values of Mcut =350, 450, 600, 
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Figure 1. 95% c.l. lower bounds on M*, as functions of m^M) for some representative values of 
Mcut (in GeV), for the four signal regions of ref. [42]. 


800, 1250, 2000, 8000 GeV. The four boxes in the figure correspond to the four different 
signal regions. The upper line in each plot, Mcut = 8 TeV, corresponds to the naive EFT 
limit, obtained without imposing any restriction on the centre-of-mass energy of the hard 
scattering.^ The limit deteriorates for decreasing Mcut because of two distinct effects. The 
first one is that the total reference cross-section g decreases, because it is restricted to a 
smaller kinematical range. This effect is unavoidable and ultimately due to the fact that 
the EFT cannot be trusted above its cutoff: trying to extrapolate the EFT above Mcut 
would be inconsistent, and this is precisely why we restrict our signal to the < Mcut 
region. The second effect is that the acceptances are also reduced, because the kinematical 
distributions of the restricted signal become softer, thus for decreasing Mcut it becomes 
increasingly difficult to pass the cuts on and on El™*®®. Being dependent on the selection, 
this effect could be mitigated by softer cuts, compatibly with the minimal El™*®® trigger 
requirement and with the fact that the SM background rapidly increases in the softer 
region. These considerations show that our signal is kinematically different from the naive 
EFT prediction: an optimized limit in all Mcut regions would require a dedicated study, 

^The naive EFT limit in SR3 differs from the ATLAS result on the D8 operator by a factor, which 
reflects the factor 2 enhancement of the cross-section for a Majorana DM particle with respect to the Dirac 
case considered in ref. [42], if the same operator is used and the normalisation in eq. (2.1) taken into account. 
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which however goes beyond the scope of the present paper and can be properly performed 
only by the experimental collaborations. 

Going back to our results in fig. 1, we notice that for large Mcut the best limits are 
obtained from the SRs with harder cuts, namely from SR2, SR3 and SR4, which all have 
comparable reach. The low-cut region SRI is instead not competitive with the other ones. 
The situation changes for low Mcut, because the cut acceptances decrease faster in the SR 
with harder cuts than in those with softer ones, and the limits start being dominated by 
the latter. For instance, when Mcut goes below 500 GeV or so, the strongest M* bound 
starts coming from SRI, while the other SRs are no longer sensitive. 

The behaviour of the limits as functions of m^M is also easily understood. When m£)M is 
lowered much below Mcut and the kinematical cuts, the cross-section becomes independent 
of rriDM and the limit saturates. The limit deteriorates as ttidm increases, because the 
latter starts having a negative impact on the energy budget of the reaction. The limit 
eventually disappears above a certain threshold, which corresponds to the region where the 
DM particle is too heavy to be produced with a centre-of-mass energy below Mcut- The 
minimal centre-of-mass energy is given by 




= + 



(2.7) 


where is the common jet and FI™'®® cut of each SR, out of which the mass threshold is 
then found to be ® 


max _ 
— 


Mcut 

2 



( 2 . 8 ) 


We thus see once again that soft SRs are favoured for low Mcut, not only because they 
produce better M* limits, but also because they have an extended reach in the DM mass.® 
The combined limits from all four SRs, obtained as the intersection of the allowed 
regions as described above or equivalently by taking the strongest M* bound at each point, 
are displayed in the left panel of fig. 2. The main conclusion we can draw is that the na'ive 
EFT limit is fairly accurate when Mcut is significantly above 1 TeV, while it considerably 
overestimates the actual exclusion for lower values of Mcut- As an equivalent way to express 
the same information, the right panel of fig. 2 shows the limit on M* as a function of 
Mcut foi' some fixed representative values of 0, 100, 250, 500, 750, 1000 GeV. This 

representation is perhaps more convenient, as the dependence on ttidm is rather smooth, 
and significant only in a limited range. Furthermore, it gives an idea of the search reach in 
the low Mcut region. For reference, the dashed line on the left-hand panel of fig. 2 shows 
the constraint from the relic density (under-abundant below the line and over-abundant 
above it), computed with the approximate analytical formulae for the EFT collected in 
Appendix G. 


®The threshold effectively occurs for lower values of tudm when Mcut gets close to the LHC threshold of 
8 TeV, because of the rapid large-x decrease of the parton distribution functions. 

^Formally, low improves the mass reach for any value of Mcut- However, at large Mcut the threshold 
has a very poor sensitivity to the actual value of pip and all SRs have practically the same reach. 
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Figure 2. Our combination of the lower bounds on M*. Left: As a function of moM, for the same 
representative values of Mcut as in fig. 1. The dashed grey line is the relic density constraint. Right: 
As a function of M^ut, for some representative values of moM (in GeV). 


The plots described above summarise the experimental situation in a simple and concise 
way, however they do not tell us how much of the theoretically allowed parameter space 
has been actually tested and how much is still unexplored. Namely, it is hard to establish 
a priori the “reasonable” Mcut values, and whether the corresponding M* limit should be 
regarded as a “strong” or a “weak” one. We can do better if we remember that Mcut and 
M* are actually connected by eq. (1.3). Clearly, we do not know what g^, is, but we do 
have some control on its value. We definitely know that it must be 5 * < dvr, since taking 
it larger would make the EFT non-perturbative below the cutoff. This implies an upper 
bound on Mcut for any given M*. In principle, there is no lower bound on g^, it could be 
arbitrarily small pushing Mcut to smaller and smaller values. However, in a WIMP-like 
scenario we definitely expect gf* ~ ~ Ij to implement the WIMP miracle recalled in the 

introduction. Values of S'* of order unity, and not radically smaller than that, should thus 
be considered as plausible benchmarks. 

The exclusion limits at fixed g^,, in the (ut-dM) Mcut) plane, are shown by the coloured 
solid lines in fig. 3, for the representative values s'* = 1.8, 2,4, 6 , dvr. The black solid line is 
the limit one would obtain in the na'ive EFT. We stress that closed excluded regions are 
obtained in this case, a fact that can be easily understood in the following terms. For a 
given moM) it is obvious that the limit must disappear at sufficiently large M*, because 
the signal cross-section rapidly decreases for increasing M*. However, the limit must also 
disappear for too low M*, because at fixed ( 7 * lowering M* means lowering Mcut = 
which deteriorates and eventually kills the signal and the acceptances. There also exist 
values of ruDM where these two competing effects do not allow to obtain an exclusion for 
any value of M*, which is why the curves close on the right. As a consequence, there are 
values of ( 7 * for which no limit on M* can be set, not even for ttidm = 0. Our finding 
is quantitatively impressive: with the experimental results available so far, a satisfactory 
exploration of the parameter space has been possible only for 5 * above 4 or 6 : the reference 
value S'* = 1 is not excluded, and the smallest coupling we are sensitive to is g^, ~ 1 . 8 . 
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Figure 3. The solid lines enclose the excluded regions in the plane (wDMjAf*), for some repre¬ 
sentative values of g*, combining the four signal regions of ref. [42]. The black line is the limit 
one would obtain with the naive EFT. The grey triangle is theoretically forbidden because of the 
self-consistency requirement M* > 2mDM/<?*, for g* = 47r. The dashed lines show, with the same 
colour code as for the solid lines, how the grey triangle expands for smaller values of g^. 


Making progress in this direction would require more energy and integrated luminosity 
at the LHC, as expected in the forthcoming runs, but also improving the sensitivity to 
the small Mcut region as explained above. Indeed, the lower exclusion limits, in the low 
rriDM region, are predicted by eq. (2.7) to occur near ~ where we 

take the lowest possible value for corresponding to 120 GeV for SRI of [42]. This 
shows once again the importance of keeping the first signal region at the lowest and 
j^miss yalues compatible with the trigger and background conditions. As a last comment, 
we remind the reader that not all the points in fig. 3 are theoretically allowed within the 
EFT framework. We are working here under the assumption of heavy-mediator DM, which 
means, as explained in the introduction, that moM should be well below Mcut, or at least 
™DM < Afcut/2, because otherwise there is no hope for the DM being produced within the 
range of validity of the EFT. This leads to the constraint M* = Mcut/9* > 9*- For 

g'* = dvr this produces the grey theoretically forbidden region in fig. 3. For < dvr the 
boundary of the grey triangle moves as indicated by the dashed lines, with specified by 
the same colour code as for the solid lines. However, eq. (2.7) guarantees that (in contrast 
with what we would obtain in the naive EFT), the experimentally excluded region can at 
most approach the theoretically excluded one. Indeed, the closeness of the solid lines to 
the corresponding dashed lines gives a measure of how much the available EFT parameter 
space has been explored for the different values of 5 *. 
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3 Simplified model reinterpretation 

In the previous section we consistently derived from experimental data universal bounds 
on the EFT dehned by the operator (2.2), as functions of the three relevant mass parame¬ 
ters (M*, moM) -^cut)- We now show how such bounds can be re-interpreted in any specihc 
microscopic model underlying the chosen effective interaction. Since it collects only the con¬ 
tribution to the (positive-dehnite) signal cross-section coming from the kinematical region 
Ecra < Mcut, where by dehnition the EFT is reliable, and it sets to zero the contribution 
corresponding to Ecm > Mcut, our prescription for using consistently the EFT leads to 
underestimating the signal cross-section. We then expect our bounds to be systematically 
more conservative than those obtained by the direct comparison of a specihc microscopic 
model with the experimental data. The aim of the present section is to perform a quan¬ 
titative comparison of the limits derived with the two methods and to comment on the 
interpretation and practical consequences of any signihcant difference in the results. 

We consider two illustrative simplihed models, characterized by quite different dynamics 
at the mediator scale, but nevertheless giving rise to the same leading effective operator (2.2) 
in the low-energy EFT. In Model A, DM annihilation into quark-antiquark pairs and the 
inverse process occur via the s-channel exchange of a spin-1 Z' boson of mass mz', coupled 
to the axial-vector currents of quarks and DM with strengths gq and gx, respectively. Very 
similar simplihed models were discussed in refs. [54-59]. In Model B, the same processes 
occur via the t/u-channel exchange of color-triplet scalars of mass fh, with the same gauge 
quantum numbers as the squarks q of supersymmetric extensions of the SM, but with a 
universal Yukawa coupling of strength goM to quarks and DM. Very similar simplihed 
models were discussed in refs. [60-65]. We have collected some useful details on the two 
models in Appendices A and B, respectively. 

Before comparing the interpretation of the experimental results in the EFT and in the 
two simplihed models, we display in hg. 4 the tree-level Feynman diagrams contributing to 
the three hard partonic processes associated with the scattering pp —?■ jet -|- MET: 

(I): q{pi) + q{p2)X{p 3 )+X{p 4 )+g{k)-, (3.1) 

(II): q{pi) + g{p2) ^ X{p3) + X{p4)+q{k); (3.2) 

(III): q{pi) + g{p2) ^ X{p3) + X{p4)+ q{k) . (3.3) 

The symbols in brackets label the four-momenta of the corresponding particles. Process I 
is described by diagram A1 in Model A, by diagrams B1 and B4 in Model B. In the case 
of diagrams A1 and Bl, it is understood that we should add the corresponding diagrams 
with the gluon radiated from the antiquark rather than from the quark line. Process II is 
described by diagrams A2 and A3 in Model A, and by diagrams B2, B3 and B5 in Model 
B, plus those obtained by exchanging the momenta pa and p 4 of the Majorana DM fermion 
X. Process III is described by the same diagrams of process II, with the prescription that 
all the arrows on the quark and squark lines should be reversed. 

The limits from our consistent EFT analysis and directly from the simplihed models 
are obtained as follows. In the EFT, we compute the EFT parameters in each simplihed 
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Figure 4. Feynman diagrams describing the jet + DM signal at hadron colliders for models 

A {Z' mediator) and B {q mediator) considered in the text. 


model and we just apply the constraints derived in the previous section. The scale M* of 
the effective operator (2.2) is given by 

M* = (Model A), M* = — (Model B). (3.4) 

yjQq 9X 9 dm 

The cutoff scale Mcnt, at which the EFT description loses its validity, is identified with the 
mediator mass Mmed) he. with mz' in Model A and with rh in Model B. Then, after this 
identification, the effective coupling g^, is: 

9* = y/9q9x (Model A), (Model B) . (3.5) 

To extract limits directly in the simplified models, we recast the ATLAS mono-jet analysis 
of ref. [42] as in section 2.1, with the only difference that now the signal cross-section is 
computed in the complete simplified model, i.e. with the diagrams in fig. 4 and with no 
Mcut restriction, for any value of M^^ed and of rriDM- For each point of the simplified model 
parameter space, the expected signal rate is computed in each SR and the corresponding 
exclusion limits are applied. 
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Model A: 95% CL limit on M, 


Model A: 95% CL limit on M, 
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Figure 5. 95% CL limit on M* for Model A, as a function of mz', for mx = 50 GeV (left) and 
mx = 250 GeV (right). The horizontal dashed line corresponds to the limit obtained in the naive 
EFT. The blue line gives the limit consistently extracted in the EFT with Mcut = mz'- All the 
other lines refer to the full model, and assume either Fz' = mz>/(87r) (solid) or Fz' = mz'/'i 
(dashed). The purple lines show the limits obtained in the full model. The red lines corresponds to 
the resonant production of the mediator. The orange lines correspond to the correct relic abundance 
for a thermal freeze-out, computed according to the formulae for Model A reported in appendix C. 
From the top left to the bottom right, the increasingly dark grey shaded areas correspond to 
Fz'Imz' > l/(87r), 1/3,1 and to g* > Att. 


For Model A, the result in the full model is illustrated by the purple lines in fig. 5, 
as an exclusion limit on M* as a function of Mmed = for two representative values 

of moM = and for two postulated values of the (width/mass) ratio of the mediator: 
Tz'lmz' = I/Stt (solid) and F^z/m^/ = 1/3 (dashed). We will see below that using the 
(niz’, M*) plane to represent the result suffers from an important limitation. Furthermore, 
M* is not a natural variable for the simplified model, where it is a derived quantity rather 
than a fundamental parameter. In this context, other ways of representing the limits could 
be more effective. The choice of the (m^/jM*) plane is however convenient for comparing 
these results with the EFT limits and with other studies of Model A, such as those in 
refs. [16, 28, 43]. In the figure, our consistent EFT limits, as reinterpreted in Model A, 
are represented by blue solid lines, while the black dashed horizontal line shows the naive 
EFT limit, formally obtained by sending Mcut to infinity for fixed M*. For reference, the 
orange lines correspond to the correct relic abundance for a thermal freeze out, computed 
here with the approximate analytical formulae for Model A reported in appendix C. 

First, we can visually check that our consistent EFT limits are actually correct model- 
independent constraints, as they lie systematically below those obtained by working directly 
with the simplified model. Notice that this is not true for the naive EFT limits, which 
overestimate the exclusion for very low mediator mass. Second, we observe that the limits 
obtained directly in Model A are slightly stronger that the EFT ones, and that this effect is 
considerably amplified for a moderately light mediator in the case of the smaller Tz'/mz' 
ratio. The reason for this behaviour is that the simplified model cross-section can get 
significantly enhanced with respect to the EFT one, leading to a stronger bound, only 
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thanks to the resonant production of the mediator, which can only take place if the latter is 
light enough. Furthermore, the resonant enhancement is of order tt mz '^-nd this is why 
it is more pronounced for a narrow mediator. These considerations are made quantitative 
by the solid and dashed red lines in fig. 5, with the same conventions as before. These 
lines represent the limits on the simplified model obtained by computing the signal rate 
restricting the invariant mass of the Z' propagators within two widths from its pole mass. 
The fact that the red lines are so close to the purple lines representing the “true” limit, when 
they are both significantly above the blue line, confirms that the resonant production is 
what drives the enhancement. It also suggests that in this kinematical region DM searches 
in the simplified model should be actually regarded as mediator searches, and the results 
reported as limits on cr{pp —)> Z')x BR(Z' —)■ XX). Also, Z' resonant production followed 
by the decay into quark-antiquark pairs, leading to a peak in the di-jet invariant mass 
distribution, may be a complementary signal to be looked for [66-69], with or without the 
extra jet: in such a case, we would obtain a limit on a{pp —)• Z')x BR(Z' —)• qq). We will 
comment further on this in the conclusions. 

We now turn to the aforementioned limitation of the plane, which was al¬ 

ready noticed in refs. [28, 43, 52], but we find important to emphasise. Model A has four 
parameters: ut-dM) W'Zq 9q, 9x- In fig- 5, the DM mass is set to a fixed value and each 
point of the plane uniquely determines mz ' and M*. Then also the product gq gx is fixed 
by the left-hand side of eq. (3.4), namely 

m|, 

9q9x — 

Only one combination of the two couplings is left free at this point, and it might seem a 
good idea to fix it point-by-point to fit the values of T z' jmz' that were assumed in drawing 
the purple lines in the figure. However, we must take into account that, for fixed gqgx^ the 
accessible values of Tz'/mz> are bounded from below: 

= (^91 +I39x> 9q9X'/^ , (3.7) 

mz' 

where a and /3 are suitably defined coefficients (see appendix A) that do not depend on 
gq and gx, and have only a mild dependence on the spectrum through phase space. This 
means that the (m^/jM*) plane is divided into regions, whose boundaries are curves (or, 
approximately, straight lines), where Fz'/mz' is always larger than a certain value. Some 
representative regions are displayed as grey shaded areas in fig. 5: from the top left to 
the bottom right, they correspond to Vz'jmz' > l/( 87 r), 1/3, 1 . The fourth and darkest 
region at the bottom right corresponds to g^ = ^/gqgx > dvr, where neither the EFT 
nor the simplified model admit a consistent perturbative description. In the neighbouring 
region where Tz'/mz' > 1, the EFT can still be consistently used, but the same does not 
apply to the chosen underlying simplified model: the fact that Tz'lmz' > 1 is telling us 
that in such strong coupling regime the simple mediator interpretation of the origin of the 
effective interaction breaks down. Even in the perturbative regime, the direct simplified 
model lines are obtained by assuming a given F^z/m^', thus they become inconsistent on 


(3.6) 
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the right of the boundary of the corresponding Vz'Irnz' region, because they cannot be 
associated to any physical point of the simplihed model parameter space. On the left plot, 
for instance, we should have stopped drawing the purple and red solid lines corresponding 
to Vz' I'mz' = l/(87r) where they cross the boundary between the white and the very light 
grey region, at mz’ ~ 600 GeV. Similarly, we should have stopped the purple and red 
dashed lines, corresponding to = 1/3, where they cross the boundary of the two 

light grey regions, at mz' ~ 1.1 TeV. The only justihcation for keeping them is that the 
limits on the width are theoretical constraints, while the actual location of the curves is the 
result of the experimental analysis, which might improve its sensitivity in the future. When 
this will happen the exclusion curves will move up and will exit more and more out of the 
inconsistent regions. As far as current data are concerned, however, this observation shows 
that the DM limits are actually rather poor, especially in the region of narrow mediator 
width, which corresponds to a weakly-interacting particle. But after all, this is exactly 
what we concluded from our exploration of the EFT parameter space: ‘small’ effective 
couplings of order one are still unconstrained. Here we have just verihed that the simplihed 
model can not help us much in this respect. 


Model B: 95% CL limit on M, 


Model B: 95% CL limit on M, 
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Figure 6. The same as in fig. 5, but for model B. The only difference is that, from top left to bottom 
right, the two diagonal lines correspond to Vqjrh = I/Stt, 1/3, and the grey areas to Tqjrh > 1 and 
to 5 * = gou!^ > dtr. 


Very similar considerations apply to Model B, whose bounds are depicted in fig. 6. Also 
in this case the enhancement of the limit obtained directly in the simplified model is mostly 
due to the resonant production of the mediator, which can occur even in the so-called ‘t- 
channel mediator’ case if an extra jet is emitted in the final state. This process corresponds 
(see diagrams B.3 and B.5 in fig. 4) to an associated DM-g production followed by the q 
decay into DM plus jet. Once again, in the region of Model B where the squarks are light 
enough, experiments should extend their selection criteria and look more generally for n = 
1 ,2 ,... jets plus to include the possibility of resonant squark production, both singly 

and in pairs. Of course, part of this is already being done in the context of standard squark 
searches within simplified supersymmetric models [70, 71]. An experimental analysis along 
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the above lines has been recently suggested by ref. [72] within a simplified supersymmetric 
model with a very light gravitino: something similar could be devised also for Model B 
and similar simplified models for DM with ‘t-channel mediators’. A second point worth 
stressing for Model B is that the issue with the (m, M*) plane is even more severe than 
in Model A, because the model has only three parameters, therefore after fixing moMi 
and M* the (width/mass) ratio of the mediator is fixed. In this case, fig. 6 shows two lines 
corresponding to Tqjfh = l/Svr, 1/3, a grey area where where Bg > m, and a dark grey 
area where > dvr. The only physical points of the four exclusion curves derived 

in Model B (purple and red, solid and dashed) are those at the intersection with the lines 
corresponding to the assumed value of Tg/m, marked as full purple dots. 

Other approaches 

We are not the first to address the issues related with the naive use of the EFT for DM in 
kinematical regimes extending beyond its range of validity: as already mentioned, they have 
been studied at length in the literature [25-33]. In particular, refs. [27, 29, 30] proposed a 
criterion (recently adopted in refs. [47, 52]) to estimate how sensitive the naive limits on 
M* are to the unsafe region of the EFT and how much they deteriorate if the latter region 
is excluded from the analysis. Below we discuss two aspects of our approach in a way that 
can be helpful for the comparison with this previous literature. 

The first point to be discussed concerns the choice of the kinematical variable to be 
used for discriminating the safe EFT region from the unsafe one. From the EFT viewpoint, 
the natural variable is clearly the hard scale of the process, E^ra'- this was our choice for 
the present paper. However, within specific ‘mediator’ models, or more precisely classes 
of models, another possible choice is the variable Qn, as proposed in refs. [27, 29, 30]. 
Qtr = +\/|Qtrl is defined as the maximal virtuality of the mediator propagator, computed 
over the Feynman diagrams contributing to the partonic DM production process under 
study. Since Qtr < Ecm, using Qtr to define the safe kinematical region of the EFT means 
gaining signal cross-section, thus obtaining a stronger and still reliable limit. Notice that, 
since the definition of Qtr depends on whether the mediator propagates in the s or in the 
t channel in the two-body annihilation qq •(-)> XX, Qtr is not suited for setting a model- 
independent limit. However, one might still consider the two possibilities in turn and set 
separate limits for the two cases of s- and f-channel mediation. While this clearly does 
not exhaust all possibilities it might be still worth doing if it considerably enhances the 
reach. 

To explore the exclusion reach of this method and compare it with ours, we start 
by recalling the (trivial) expressions for Qtr in Models A (s-channel) and B (f-channel), 
corresponding to the diagrams in fig. 4 and the conventions in eqs. (3.1)-(3.3). In Model A, 
for both process I and process H (the kinematics of HI is identical to that of H, so it does 
not need a separate discussion), Qtr is just the invariant mass of the DM pair 

Qtr = iP3 + Pif = {pi +P2- kf (AI, AH) . (3.8) 

^°The effective interaction might well be generated by the combined exchange of s- and t-channel medi¬ 
ators, or by radiative effects not falling in any of these two categories. 
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In model B, instead, we have to consider process I and II,III separately. In the case of 
process I, Qtr reads 

Ql = max {{pi-k- = (ps - P 2 )^ {pi - Pif = {Ps - P 2 + kf , 

{pi - k - psf = {pi - P2f , {pi - ps)^ = {pi - P2 + kf} (BI), (3.9) 

while for process II,III we have 

Qtj. = max {(pi -psf = (p 4 - p 2 + kf , (pa + kf = (pi + P 2 - P4)^ , 

{Pi-Pif = {P3-P2 + kf , {pi + kf = {pi+P2-p3f} (BII). (3.10) 

Notice that the subprocesses are quantum-mechanically distinguishable and therefore it 
makes sense to adopt a different definition of Qtr for each of them. 




Figure 7. Limits on M* as functions of M^md obtained for Models A and B with three different 
methods. The purple lines are derived in the full models, assuming two representative values of the 
ratio Tnied/Afnied: l/( 87 r) (solid) and 1/3 (dashed). The solid blue line is derived in the EFT with 
our method as described in the text. The solid green line is derived in the EFT by imposing the 
condition on Qtr proposed in refs. [27, 29, 30]. Upper plots: Model A. Lower plots: Model B. 

^^Notice that, if in Model B we had assumed a Dirac DM particle, only half of the conditions in eqs. (3.9) 
and (3.10) should have been imposed. Therefore, the model dependence of this strategy depends on the 
assumptions made both on the mediator (s-channel or t-channel) and on the nature of the DM particle 
(Dirac or Majorana fermion, complex or real scalar). 


- 18 - 




















The result of the comparison is displayed in fig. 7, where we show the limits on M* 
as functions of Mmed) obtained for Models A (upper plots) and B (lower plots) with three 
different methods. The purple and blue lines represent the full model and our approach to 
the EFT, respectively, namely the same curves as in figs. 5 and 6. The green line is also 
derived in the EFT, but with the cut Qtr < instead of Elcm < Mcut- In the limit 

of heavy mediators, all the lines coincide as expected. The differences are in the region of 
relatively light mediators, where the EFT limit obtained with Qtr has, as expected, a better 
reach in M* than our method. However, in our view the improvement is not sufficiently 
significant, especially when compared with that obtainable in the full simplified model, 
to motivate the use of Qtr rather than Ecm- Our recommendation is thus to stick to the 
simple and model-independent version of our method, possibly trying to extend the reach 
by the direct search of the mediator which, as described in the previous section, is the sole 
responsible of the improved reach of the simplified model. 

The second aspect to be clarified is that the consistent EFT limits in the [rriDM, M^) 
plane, at fixed g^, cannot be inferred from those obtained in the na'ive EFT by just per¬ 
forming a rescaling of M*. It is indeed clear that such a rescaling cannot lead to closed 
exclusion curves such as those we obtained in fig. 3. One might be tempted to consider a 
rescaling here because the EFT cross-section is proportional to l/M^, so that the reduction 
of the cross-section caused by the kinematical cut might be reabsorbed into an effective M*. 
Namely, one might consider defining the ratio 




a{M^) 


Qtr'^^cut — 

a{M^) 


(3.11) 


where a denotes the signal cross-section computed in the na'ive EFT for a given signal 
region. At fixed and tudM) 72 is a function of M*, which tends to one for sufficiently 
high M* and to zero for sufficiently low M*, because of the effect of the kinematical cut 
illustrated in eq. (2.7). Given that R measures the reduction of the cross-section with 
respect to the na'ive EFT, one might think of getting the limit on M* at each tudM; call 
it M*, starting from the limit obtained in the na'ive EFT, call it and solving the 

implicit equation 


M* 


72 M*, moM 


i ^EFT 


(3.12) 


The effective operator scale M* obtained in this way is the one that reproduces, in the EFT 
with the cut on Qtr, the same signal cross-section that was needed for setting the bound at 
na'ive EFT. Namely, eq. (3.12) is equivalent to 


cj(MQ 


Qtr<Mcut=g*M* 


aiM, 


EFTt 


(3.13) 


where we have exploited the fact that in the na'ive EFT a{M^) simply scales as 1/M^. 

The above method for obtaining M* is more elaborate than directly comparing the 
experimental limit on the cross-section with the prediction of the kinematically restricted 

^^Using Qtr or Ecm makes no difference for the point we want to make here. 


19 - 








EFT, as we suggested in section 2.1. Furthermore, the rescaling method might obscure the 
fact that eq. (3.13), or equivalently eq. (3.12), has either zero (which means that no limit can 
be set) or two solutions for M*, but it never has only one. The behaviour of the restricted 
EFT cross-section, compared with the naive EFT, is pictorially represented in fig. 8 . The 



Figure 8 . A pictorial representation of how the bounds on depend on the prescription for 
computing the signal in the EFT. The signal cross-section is displayed as a function of M*, for fixed 
g* and moM- The black and the blue lines correspond to the naive EFT and to our consistent pre¬ 
scription, respectively. The horizontal purple line represents the experimental limit. The resulting 
excluded interval for M* is reported near the horizontal axis for the two prescriptions. 

cross-section vanishes before approaching M* = 0, because of the cut Qtr < g*M^. There¬ 
fore there are two values of M* for which the cross-section equals the experimental limit, 
which means that the excluded region has one upper but also one lower limit in M*, dif¬ 
ferently from the naive EFT as depicted in the figure. Therefore, the true limit cannot be 
set by just rescaling the naive EFT exclusion curve. Notice also that certain strategies to 
solve eq. (3.12), such as applying an iterative procedure starting from M* = might 

obscure the existence of the lower bound, as they systematically converge to the upper 
one. The quantitative impact on the excluded regions in the (ut-dm, ^*) plane, for different 
values of g^,, was already displayed in fig. 3 for our kinematical requirement Ecm < g*M^. 

4 Conclusions 

We described a simple strategy to set robust and mo del-independent limits on heavy- 
mediator DM at the LHC. Our method is based on the generic form of the operators in the 
EFT containing only the DM and the SM particles, with no assumptions on the underlying 
dynamics. However, it also takes into account the presence of a cutoff scale Mcut above 
which the EFT loses its validity. M^ut must be regarded as one of the free parameters of 
the EFT, on the same footing as the DM mass ttidm and the effective interaction scale M*. 
We have to do so if we aim at a comprehensive exploration of the whole range of theoretical 
possibilities. The parameter M^nt can be traded for the typical coupling strength at the 
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mediator scale. As explained in the paper, can be defined in the EFT alone, and further 
characterised for any assumed underlying model. 

We applied our method explicitly to the ATLAS mono-jet search of ref. [42], obtaining 
the exclusion contours in the (ttidM) Af*) plane shown in fig. 3, for fixed representative values 
of ( 7 *. We believe that this kind of plots illustrates the current experimental situation in 
an accurate and comprehensive way, providing a fair assessment of the LHC sensitivity 
to heavy-mediator DM. At the moment, we are only sensitive to large values of g^, while 
the region s'* ~ 1, which is arguably the most natural one for WIMP DM, is still largely 
untested. Making progress requires higher energy and luminosity, but also an optimisation 
of the experimental search strategies. As pointed out in section 2, our signal is kinematically 
different from that of the naiVe EFT, in particular it is characterised by softer and 
distributions. The reach of the searches would then benefit from a sensitivity improvement 
in the soft region. 

In section 3 we considered two different simplified models, both giving rise to the same 
effective operator considered in section 2. We compared our EFT limits, reinterpreted in the 
two models, with those obtained from a dedicated comparison of the experimental bounds 
with the prediction of the two models, reaching two main conclusions. First, the limits 
set within the simplified models can be considerably stronger than the EFT ones, but only 
because of the resonant production of the mediator, which enhances the simplified model 
cross-section. Therefore, a DM search performed within a simplified model (in the only 
interesting region where the limit is potentially stronger than in the EFT) is actually not 
a search for DM, but a search for the mediator, and as such it should be interpreted. The 
canonical fixBR limit as a function of the mediator mass appears to be the best option for 
presenting the experimental results. The second conclusion is that the current experimental 
sensitivity is still rather poor, even when working within a simplified model. In particular 
the region of weak coupling, i.e. narrow mediator, is mostly unexplored, in accordance with 
what we found in our EFT analysis. We finally discussed two aspects of our approach, to 
facilitate the comparison with the recent literature. We found that the usage of the variable 
Qtr hr place of to define the safe kinematical region of the EFT does not improve the 
sensitivity significantly enough to pay back for the increased model-dependence. We also 
remarked that just rescaling the naive EFT limit does not account for the impossibility, 
within mono-jet searches, of excluding arbitrarily low values of M* at fixed rriDM and g^,. 

In summary, we have found that the LHC sensitivity to the heavy-mediator DM hy¬ 
pothesis is still limited and wide regions of the parameter space still wait to be explored. On 
the experimental side, improving the analysis in the soft region would be of great help. On 
the phenomenological side, more comprehensive methods should be elaborated to cover each 
different region of the parameter space with the most suitable strategy. Non-resonant DM 
signals are well described by the EFT which, as outlined in the present paper, when con¬ 
sistently used provides a robust model-independent way to approach the problem. Within 
specific models, this needs to be supplemented by resonant mediator searches, which how¬ 
ever should be performed by exploiting fully the predictive power of the assumed mediator 
dynamics. This means taking into account all the mediator production mechanisms (single 
and/or pair) and all its possible decay modes, including the one to visible objects which 
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might give complementary informations. These aspects are left to future work. 
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A Model A: axial-vector mediator 


We collect here some details on the first of the two simplified models considered in the text. 
Model A. Previous discussions of very similar models can be found in refs. [54-59]. The 
mediator is a neutral vector boson Z', singlet under the SM gauge group, with mass mz', a 
universal axial coupling Qq to quarks, no renormalisable couplings to leptons, and an axial 
coupling gx to the Majorana DM fermion X of mass mx- Since the model is introduced for 
purely illustrative purposes, without making reference to an underlying more fundamental 
theory, we introduce an explicit Z' mass term and we neglect Z-Z' mixing, as well as 
anomalies and their cancellation mechanisms. 

The model Lagrangian is 


Ca — Fsm + dlx + dlz' + 1 



(AT) 

(A.2) 

(A.3) 


where Csm is the SM Lagrangian, Cx is the free Lagrangian for X in (2.1), Z'^^^ = — 

dyZ'^, and the sum in (A.3) runs over all quark flavours {q = u, d, c, s, t, b). The model has 
four parameters, 

gq, gx, = Mined, mx ^ mr,M, (A.4) 


which can all be taken to be real and positive (in principle, gq and gx could have either 
sign, but this is not relevant for the present study). 

Notice that the choice of a purely axial interaction, universal for all quark flavours, is 
crucial to generate the effective interaction (2.2) from (A.3) in the low-energy limit. At 
leading order in E/mz' ^ 1, the approximate solution of the Z' equations of motion is 


= -- 




- 

Z' 1 
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which substituted in (A.l) gives 





(A.5) 

q q 

(A.6) 

771 ry / 

(A.7) 


The effective interaction term (A.7) between the SM quarks and the DM field reproduces 
the one in (2.2) as long as 

M, = . (A.8) 

^/9q9X 

Notice also that integrating out the heavy Z' generates two additional four-fermion op¬ 
erators, (A.5) and (A.6). However, (A.5) is subject only to very mild constraints from 
the limits on DM self-interactions. The four-quark operator (A.6) can be probed by the 
searches for contact interactions [73, 74], but can be parametrically suppressed by choosing 
9X > 9q for fixed g^. 

At tree-level, and including only two-body decays, the total decay width of the Z' is 

.2\3/21 
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rrir 
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mr 


(A.9) 


with the obvious modifications if some of the final states are not kinematically accessible. 

B Model B: coloured scalar mediators 

We collect here some details on the second of the two simplified models considered in the 
text. Model B. Previous discussions of very similar models can be found in refs. [60-65]. 
In Model B, the interactions between the SM quarks and the DM particle X are mediated 
by three families of degenerate complex scalars of mass m, with the same gauge quantum 
numbers of the corresponding left- and right-handed quarks. Since they are identical to the 
squarks of supersymmetric extensions of the SM, we denote them with the same symbols, 
(uiL, diL,Uiii, dm), where i = 1, 2, 3 are family indices. Similarly, the Majorana fermion X 
mimicks, although in the simplified fashion specified by its interactions below, the lightest 
neutralino of supersymmetric models. 

The model Lagrangian reads 


Cb = CsM + ^X + ^q +^fnt ^ (B.l) 

{d>"uiL)\d^UiL) + (d^diL)^ {d^dii) + {d^UiR}^ {df^UiR) + id^diR)\d^diR) 


3 

a = E 

2=1 


- {uiL^iL + d^^diL + ul^UiR + d^j^diR^ 


+ 


^int 9 dm 


(uiL UiL + diL diL + UiR UiR + diR dm'^ X + h.c. 
1=1 ^ 


(B.2) 

(B.3) 
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where Csm and Lx are the same as in Model A, and the dots in (B.2) denote the squark 
gauge interactions, generated by promoting ordinary derivatives to SM covariant derivatives. 
Notice that the mass degeneracy and the universality of the Yukawa couplings among 
quarks, squarks and DM evade the typical problems of supersymmetric models with flavour¬ 
changing neutral currents. The model has three parameters. 


Qdm , m — Mmed , nrx = rriDM , 


(B.4) 


which can all be taken to be real and positive {goM can be complex, but it can be chosen 
to be real and positive by absorbing its phase into a redehnition of the squark helds). 

As for Model A, we can derive the EFT by solving the classical equations of motion 
for the squarks in the low-energy limit if <C m: 


fi'DM ^ 

UiL = - A UiL , UiR = 


Substituting into Lb yields 


rB _ 
‘-EFT — 


ffpM 
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[{XuiL){uILX) + {XuiR){lHkX) + {XdiL){d~LX) + 
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gpM 
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■ 3 

-i=l 


(B.6) 


where for the second equality we have used the Fierz identities and the fact that when A is 
a Majorana spinor Ay^A = 0. The effective interaction term (B.6) between the SM quarks 
and the DM particle reproduces the one in (2.2) as long as 




2 m 
Qdm 


(B.7) 


At tree-level, and assuming m > mx + Tfiq, where q is the corresponding quark, the 
decay width of the generic q is 


IGtt 


m 2 
5pm 1 


(m^ -|- 


^ I V -9 ^ ruq + mx 

nqA ^2 



(B.8) 


C Formulae for the relic density 

We collect here the approximate analytical formulae used for the calculation of the relic 
density in the EFT (hg. 2), in Model A (hg. 5) and in Model B (hg. 6), before requiring 
that it reproduces the recent precise determination by the Planck collaboration [75] (for 
our purposes, the latter can be rounded to Vlj)y[h? = 0.12 with negligible error). They can 
be straightforwardly derived from the existing literature (see e.g. [76, 77]). Up to terms of 
order 1/x/, where Xf is the value of x = mx/T at freeze-out: 

- 1-07 • 10^ (GeV)-^ , (C.l) 
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where h is the dimensionless Hubble parameter, ~ 100 is the number of relativistic 
degrees of freedom, Mp ~ 2.4 x 10^® GeV is the reduced Planck mass, mx is the mass of 
the DM particle in GeV, and 


Xf = ln(A) — - ln[ln(A)] + In 


A = 0.038—;—APp uix 


1 + 


66 1 


a ln(A)J ’ 


lQm\ 


In the EFT, introducing the dimensionless parameters aq = niq/mx, 


a = 


E 


2 /. _ 2 

TT I mJ 




(C.2) 
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where the sums run over the quark flavours whose mass is below mx- 

In the two models underlying the EFT, we introduce two additional dimensionless 
parameters, j3 = mx/M^ned and 7 = Tmed/A^med) to account for the hnite mass M^ed and 
width Tjned of the mediator. Then in Model A (Z' mediator) 



(1 - 8/32 + 16/3^) , 
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and in Model B {q mediator) 
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